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RESEARCH  OVERVIEW 


This  report  covers  the  period  from  October  1,  1985  through  March  31,  1986. 

The  research  discussed  here  is  described  in  more  detail  in  several  published 
and  unpublished  reports  cited  below. 

During  this  period,  the  installajLijm.  of  the  .ceatcal.cfiJDpittiiur  machine  for  the 
MIT  CAP  work,  known  as  CAP 'Tor  more  formally,  as  caf .mit.edup^and  all  the 
necessary  software  was  substantially  completed.  The  machine  is  a  VAX  11/785. 
It  has  a  large  cadre  of  users. 

A  rudimentary  CAP  system,  known  as  CAPE^ (Computer-Aided j^PabricationJEnviron- 
ment),  has  been  implemented.  So  far  this  system  ^ovid«  eui  electronic  per¬ 
sonal  notebook,  machine  instructions,  and  equipment  status  reports.  Four  fab 
rication  facilities  are  being  served,  at  least  in  part.  Even  though  crude, 
CAFE  has  been  in  constant  use. 

The  requirements  of  a  scheduling  program  are  being  studied.  One  surprising 
result  so  far  is  the  importance  of  setup  operations  in  scheduling  some  equip¬ 
ment. 


Two  new  process-simulation  modules  have  been  developed.  One  is  for  a  focused 
ion  beam  operation  used  for  implanting,  and  the  other  is  for  nonuniform  oxida 
tion.  Improved  models  were  inserted  into  the  simulators  SUPREM  and  MINIMOS. 


SCHEDULING 


Research  during  this  period  focused  on  three  activities:  studying  the  inte¬ 
grated-circuit  fabrication  process  at  a  systems  level,  formulating  a  mathema¬ 
tical  model  of  an  integrated-circuit  fabrication  facility,  and  developing  an 
electronic  sign-up  sheet  for  scheduling  the  laboratory. 

The  first  study  is  to  define  the  scheduling  problem.  This  involves  studying 
VLSI  technology,  the  construction  of  the  MIT  Integrated  Circuits  Laboratory, 
and  user  requirements.  It  was  discovered  that  a  few  key  elements  of  the  prob¬ 
lem  may  have  dramatic  impact  on  the  final  schedule.  A  setup  is  required  when 
there  is  an  important  change  in  type  of  operation  at  a  machine,  and  can  be  due 
to  changes  of  gases,  water  temperatures,  or  other  required  materials.  Setups 
complicate  scheduling  because  machines  that  are  undergoing  such  changeovers 
are  not  available  for  productive  operations.  There  is  an  incentive  to  put  off 
changeovers  as  long  as  possible  because  when  they  take  place  often,  production 
capacity  is  reduced.  On  the  other  hand,  if  changeovers  are  performed  too 
infrequently,  some  lots  will  be  delayed  a  long  time. 

Randomness  is  due  to  machine  failures  and  unpredictable  demand  on  the  facil¬ 
ity.  In  a  university  setting,  the  mix  of  chips  going  through  the  facility  has 
very  little  consistency  from  week  to  week.  We  will  have  to  group  the  demands 
for  scheduling  purposes,  and  use  guesses  based  on  historical  data.  The  pres¬ 
ence  of  randomness  means  that  any  schedule  must  have  some  idle  time  in  it,  so. 
that  a  perturbation  does  not  propagate  indefinitely.  It  also  means  that  there 
must  be  some  means  for  rescheduling  after  events  like  machine  failures. 

We  are  also  investigating  the  differences  between  a  university  laboratory  and 
a  commercial  production  system.  In  a  commercial  facility,  machine  failures 
are  no  less  important,  but  there  is  a  more  predictable  demand  on  the  facility. 
The  CAF/Scheduling  group  visited  DEC  in  Hudson,  Massachusetts  on  February  11, 
1986  to  observe  their  production  system. 

The  second  activity  is  to  develop  a  scheduling  model  of  the  system.  Develop¬ 
ing  such  a  model  requires  identifying  scheduling  objectives,  decomposing  the 
problem,  solving  each  of  the  subproblems,  developing  a  database  to  store  the 
required  information,  and  testing  the  solution  on  the  facility.  Scheduling  is 
complicated  by  setups  and  randomness,  and  various  sources  of  both  are  being 
studied.  The  main  focus  of  the  analytic  activity  is  the  setup  issue.  The 
approach  now  being  explored  views  setups  as  events  that  take  place  at  a  lower 
frequency  than  machine  operations  and  most  other  events.  A  slow- time-scale 
model  is  being  explored  that  accounts  for  the  time  that  the  system  spends  in 
each  setup  configuration,  as  well  as  the  time  spent  in  changing  configura¬ 
tions.  The  impact  of  setups  on  capacity  is  modeled  explicitly. 

Remaining  work  on  this  formulation  includes  modeling  the  scheduling  objective, 
and  then  translating  the  setup  frequencies  that  are  calculated  this  way  into 
actual  times  at  which  to  perform  the  changes  in  configurations. 


The  purpose  of  this  activity  is  to  develop  a  simple  Version  1  scheduler  which 
will  not  use  the  sophisticated  ideas  that  are  described  above,  but  which  can 
be  implemented  and  used  relatively  quickly.  The  existing  paper  system  is 
being  designated  as  Version  0,  and  the  first  electronic  system  as  Version  1. 
Great  consideration  is  being  given  to  ease  of  use.  It  will  be  designed  to 
emulate  the  paper  sign-up  sheet  in  the  old  facility,  with  no  additional 
features  except  for  the  fact  that  it  can  be  accessed  from  any  terminal. 
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MODULAR  PROCESS 


The  focus  of  this  project  is  the  development  of  methods  for  IC  fabrication 
process  and  device  design  in  the  context  of  a  Computer  Aided  Fabrication  (CAP) 
research  laboratory.  The  main  thrust  of  our  work  centers  on  the  development 
of  what  we  call  the  MASTIF  workstation  (MIT  Analysis  and  Simulation  Tools  for 
IC  Fabrication).  This  is  a  menu-  and  window-oriented  program  written  in  C  and 
FORTRAN  which  provides  methodology  and  uniform  software  structure  for  the 
connection  of  process  and  device-design  CAD  tools.  MASTIF  currently  includes 
a  facility  for  incremental  development  and  version  management  of  a  process 
description,  management  mechanisms  for  definition  of  physical  cross  sections 
deriving  from  the  overall  process  description,  and  an  interactive  graphics 
interface  and  data  interchange  for  process  and  device  simulators  (SUPREM-III 
and  MINIMOS).  With  MASTIF  the  user  can  effectively  develop  and  evaluate  a 
fabrication  process  via  a  single  integrated  workstation. 

In  parallel  with  the  development  of  MASTIF,  work  is  under  way  in  the  area  of 
process  and  device  modeling.  Extensive  use  of  SUPREM-III  has  revealed  several 
shortcomings  of  the  program.  Perhaps  the  most  significant  improvement  that  we 
have  implemented  in  an  experimental  fashion  is  the  inclusion  of  dynamic  clus¬ 
tering  for  arsenic  atoms  in  silicon.  Similarly,  we  have  scrutinized,  found 
inaccurate,  and  improved  significantly  the  electron-mobility  model  in  the 
device  simulator  MINIMOS. 

Two  new  process-simulation  modules  have  been  developed  during  this  period. 

The  first  is  used  for  the  generation  of  as-implanted  and  diffused-doping  pro¬ 
files  generated  by  means  of  Focused  Ion  Beam  (FIB)  implantation.  This  is  in 
support  of  a  separate  FIB  applications  research  program  also  supported  by 
DARPA. 

The  second  process-simulation  module  is  still  in  an  experimental  form.  It 
allows  the  simulation  of  nonuniform  silicon  oxidation.  The  calculated  results 
are  the  shape  of  the  grown  SiOn  and  the  stress  exerted  during  the  process  at 
the  Si/Si02  interface.  A  boundary-integral-equation  method  has  been  developed 
for  the  numerical  solution  of  the  problem.  Because  of  the  many  uncertainties 
in  physical  parameters  for  grown  Si02,  the  present  simulator  allows  the  inves¬ 
tigation  of  three  different  oxide-movement  behaviors:  elastic,  viscous,  and 
visco-elastic.  Work  is  under  way  to  establish,  by  comparison  with  experimen¬ 
tal  data,  the  range  of  validity  of  these  different  models. 


CAF  SYSTEM  STRUCTURE 


Prior  to  the  beginning  of  this  reporting  period,  an  initial,  rudimentary  CAF 
system  was  designed  and  implemented.  This  included  a  personal  notebook  and 
data  structures  for  the  upcoming  fabrication  equipment  installation.  This  vas 
implemented  on  a  borrowed  VAX  11/750  and  some  of  the  engineering  staff  started 
to  use  the  lab  notebook  and  editor.  RS232  cables  to  the  offices  and  labs  were 
installed  and  some  were  checked  and  connected  to  terminal  concentrator  ports. 
Three  concentrators  were  assembled,  connected  to  an  MIT  Ethernet,  and  were 
operational.  The  VAX  11/785  hardware  was  installed. 

Since  then,  we  have  completed  the  computer  hardware  installation,  evaluated 
and  purchased  terminals  for  both  users  and  for  use  in  the  clean  rooms, 
installed  Unix  version  4.3BSD,  integrated  local  networking  software,  adapted 
mail  software  to  our  local  needs,  ported  several  applications  programs  to  our 
CAF  computer,  instituted  a  weekly  back-up  procedure,  substantially  improved 
our  terminal  concentrators,  created  user  documentation,  interfaced  to  the  MIT 
physical  plant  computer  which  monitors  a  large  number  of  sensing  points  in  our 
facility,  initiated  some  student  projects,  initiated  plans  for  including  the 
teaching  laboratories  in  our  developing  CAF  system,  acquired  both  documenta¬ 
tion  and  software  for  the  present  Berkeley  CAF  system,  initiated  acquisition 
of  additional  main  memory  and  disk  storage,  and  expanded  our  user  base  to 
include  virtually  all  of  the  faculty,  students,  and  technical  staff  associated 
with  our  fabrication  facility. 

Our  present  computer  hardware  consists  of  a  VAX  11/785  with  8  Mbytes  of  main 
memory,  two  450  Mbyte  disks,  GCR  tape,  1600-cpi  tape,  laser  printer,  four 
phone  lines  and  modems,  and  an  Ethernet  port  to  the  MIT  network.  Additionally 
there  are  seven  11/73-based  terminal  concentrators,  each  of  which  has  32  RS232 
ports.  RS232  cables  have  been  installed  throughout  our  building,  and  numerous 
user  terminals  and  PCs  have  been  connected. 

We  have  completed  an  extensive  evaluation  of  terminals  for  use  in  the  clean 
rooms,  and  by  the  technical  staff  and  other  CAF  users.  The  information 
gleaned  during  our  selection  process  has  been  communicated  to  other  building 
occupants  for  their  consideration.  We  chose  Ergo  301  terminals  as  they  had 
the  best  compromise  considering  cost,  video  quality,  compatibility  with 
existing  software,  and  the  fact  that  they  do  not  have  fans,  which  allows  their 
use  in  the  clean  room  areas. 

We  installed  the  latest  version  of  UNIX,  Berkeley  4.3BSD,  and  have  integrated 
local  networking  software  and  have  adapted  the  mail  software  to  enable  local 
as  well  as  remote  communications.  Our  official  host  name  is  caf.mit.edu; 
locally  the  computer  is  known  simply  as  caf. 

Several  large  applications  programs  have  been  successfully  ported  to  run  on 
mit-caf,  including  MAGIC  and  PISCES.  The  success  in  porting  these  programs 
has  resulted  in  greater  use  of  the  caf  system,  so  much  that  we  will  soon  need 
increased  primary  memory  and  user  disk  space. 
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A  back-up  regimen  has  been  initiated  with  full  dumps  taken  every  month  and 
incremental  dumps  every  week.  This  is  intended  primarily  as  protection 
against  catastrophic  disk  failure. 

Experiments  with  our  terminal  concentrators  revealed  service  deficiencies  when 
all  ports  were  active.  The  cause  was  traced  to  inadequate  use  of  the  avail¬ 
able  buffer  memory  by  the  SWITCH  software.  This  has  now  been  corrected  and 
satisfactory  performance  is  achieved  with  32  active  lines  per  concentrator. 

We  are  reasonably  confident  that  we  could  implement  48  lines  per  concentrator 
but  are  not  sure  that  satisfactory  response  times  would  be  preserved  if  all  48 
lines  were  active  at  once.  We  have  more  physical  memory  than  can  now  be  used, 
probably  enough  to  service  twice  the  number  of  lines  per  concentrator.  How¬ 
ever,  it  is  doubtful  that  the  rather  substantial  software  rewrite  required  to 
use  this  memory  is  worth  the  effort. 

We  have  automated  the  procedures  for  applying  for  and  approving  new  computer 
accounts.  An  program  named  "open"  has  been  written  which  enables  prospective 
users  to  make  requests  for  accounts  and  provide  the  required  information. 

Mail  is  sent  to  an  "approver"  after  an  account  has  been  requested  and  new 
accounts,  including  appropriate  initialization  files  and  directories,  are 
automatically  created  after  approval. 

We  have  also  been  active  in  the  generation  of  documentation  for  CAP  users. 

User  documentation  now  covers  the  use  of  terminal  concentrators,  obtaining 
accounts,  obtaining  documentation,  using  UNIX  and  EMACS,  using  Kermit  with  IBM 
PCs,  and  using  the  nroff  text  formatter  with  the  -me  macro  package.  Manual 
entries  for  the  Cshell  with  command  completion  and  Mail  utility  have  also  been 
made  available  for  distribution.  We  have  initiated  a  project  to  define  laser 
printer  hardware  and  software  to  provide  both  typesetting  and  hard  copy  of 
arbitrary  graphics. 

We  have  interfaced  to  the  MIT  physical  plant  computer  which  monitors  a  large 
number  of  sensing  points  such  as  resistance  of  DI  water,  etc.,  in  our 
facility.  We  capture  all  alarms  printed  and  have  provided  a  mechanism  for 
remote  query  of  the  status  of  the  monitored  points.  Our  facilities  personnel 
can  now  interrogate  the  monitoring  computer  from  either  office  or  home  and 
efficiently  direct  corrective  action. 

We  have  initiated  a  plan  to  expand  the  availability  of  the  CAP  system  to  the 
integrated-circuit  teaching  laboratory.  This  will  provide  us  with  another 
active  test  bed  for  the  CAP  system,  and  will  also  allow  indoctrination  of 
entering  graduate  students  (our  present  practice  is  to  require  new  students  to 
complete  an  intensive,  short  version  of  the  undergraduate  teaching  laboratory 
subject) . 

We  have  advertised  to  our  students  several  projects  concerned  with  integrating 
measuring  instruments  into  the  CAP  system.  One  student  has  started  work  on  a 
project  relating  to  a  CV  measuring  instrument.  Another  project  has  been 
defined  concerning  software  for  a  scanning  thermometer  intended  for  use  in 
monitoring  photolithography  process  parameters. 
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Ve  have  acquired  both  documentation  and  software  for  the  present  Berkeley  CAP 
system.  Ve  are  in  the  process  of  installing  a  test  version  of  this  software 
and  will  study  the  integration  of  this  software  into  our  system  or  the  inte¬ 
gration  of  our  system  into  the  Berkeley  software. 

Ve  now  have  a  substantial  number  of  CAP  users.  The  system  response  time  re¬ 
mains  quite  good,  but  we  have  had  to  make  plans  to  augment  primary  memory  and 
available  disk  space.  Our  primary  future  activity  is  to  further  the  develop¬ 
ment  of  the  CAP  applications  software  and  its  use  in  the  installation  and 
operation  of  equipment  to  be  installed. 


8 


PUBLICATIONS  LIST 


T.-L.  Tung  and  D.  A.  Antoniadis,  "A  Boundary  Integrated  Equation  Approach  to 
Oxidation  Modeling,"  IEEE  Transactions  on  Electron  Devices,  vol.  ED-32,  no. 
10,  pp.  1954-1958,  October  1985. 

D.  Boning  and  D.  A.  Antoniadis,  "MASTIF  -  A  Workstation  Approach  to  Fabrica¬ 
tion  Process  Design,"  Digest  of  Technical  Papers,  IEEE  International  Confer¬ 
ence  on  CAD,  ICCAD-85,  Santa  Clara,  CA,  pp.  280-282,  November  18-21,  1985. 

C.  Lozinski  and  S.  B.  Gershwin,  "Dynamic  Production  Scheduling  in  Computer- 
Aided  Fabrication  of  Integrated  Circuits,"  to  be  presented  at  the  IEEE 
Robotics  and  Automation  Conference,  San  Francisco,  CA,  April  7-10,  1986. 


TALKS  WITHOUT  PROCEEDINGS 


P.  Penfield,  Jr.,  "Computer-Aided  Fabrication  and  Other  VLSI  Research  at  MIT 
IBM  Corporation,  Yorktown  Heights,  NY,  March  21,  1986. 

D.  Boning  and  D.  A.  Antoniadis,  "MASTIF  -  A  Workstation  Approach  to  Fabrica¬ 
tion  Process  Design,"  MIT  Fall  1985  VLSI  Research  Review,  Cambridge,  MA, 
December  16,  1985. 


1954 


IEEE  TRANSACTIONS  ON  ELECTRON  DEVICES.  VOL.  ED-3:.  NO.  10.  OCTOBER  1985 


A  Boundary  Integral  Equation  Approach  to  Oxidation 

Modeling 

THYE-LAI  TUNG,  student  member,  ieee,  and  DIMITRI  A.  ANTONIADIS,  senior  member,  ieee 


Abstract — Thermal  oxidation  of  silicon  involves  the  diffusion  of  oxi¬ 
dant  molecules  from  the  gas-oxide  interface  to  the  oxide-silicon  inter¬ 
face,  and  the  transport  of  newly  formed  oxide  away  from  the  latter. 
Under  suitable  formulations  these  two  processes  can  be  shown  to  be 
boundary-value  problems  of  harmonic  and  biharmonic  nature.  Based 
on  these  properties,  a  boundary  integral  equation  method  (BIEM)  has 
been  developed  for  modeling  oxidation.  This  method  achieves  simplicity 
and  efficiency  by  solving  a  two-dimensional  problem  using  line  integrals 
on  the  boundaries.  The  use  of  source  distributions  as  intermediary  so¬ 
lutions  facilitates  direct  calculations  of  a  wide  variety  of  boundary  pa¬ 
rameters. 


Introduction 

IN  local  oxidation  of  silicon  (LOCOS)  process,  the  shape 
of  the  oxide  undergoes  drastic  transformation  because 
of  nonuniform  oxide  growth  rate.  A  numerical  method  that 
can  model  such  change  with  ease  is  desirable  for  compre¬ 
hensive  IC  fabrication  simulation.  Methods  based  on  in¬ 
tegral  equations  have  previously  been  explored  for  mod¬ 
eling  oxidation  [1],  [2].  Compared  to  conventional 
methods  such  as  finite  difference  and  finite  element,  they 
have  achieved  some  advantages,  for  instance,  larger  tol¬ 
erance  on  the  deformations  of  the  mesh;  nevertheless  they 
still  require  calculations  within  the  oxide  bulk. 

This  paper  reports  a  Boundary  Integral  Equation 
Method  (BIEM)  for  modeling  thermal  oxidation  in  two 
dimensions.  Our  method  relies  on  Green’s  theorem  for 
harmonic  and  biharmonic  problems,  and,  unlike  all  pre¬ 
vious  approaches,  obtains  solutions  solely  from  line  inte¬ 
grals  along  the  oxide  boundaries,  achieving  reduction  in 
dimensionality.  Because  it  does  not  require  a  mesh  to  sub¬ 
divide  the  oxide  bulk,  it  is  simple  and  efficient.  We  only 
have  to  deal  with  the  boundaries,  precisely  what  is  needed 
in  many  applications. 

Physical  Models 

The  two  constituent  processes  of  thermal  oxidation  are 
the  diffusion  of  oxidant  and  the  motion  of  oxide.  In  the 
oxidant  transport  process,  oxidant  molecules  diffuse  from 
the  ambient,  through  the  oxide  bulk,  to  the  oxide-silicon 
interface  to  react  with  silicon.  Due  to  its  volume  expan¬ 
sion.  newly  created  oxide  partly  fills  the  void  left  by  con¬ 
sumed  silicon,  and  partly  displaces  and  pushes  existing 
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oxide  towards  the  surface.  At  sufficiently  high  tempera¬ 
tures,  oxide  behaves  like  a  viscous  fluid,  thus  its  motion 
can  be  studied  as  a  fluid  flow  problem. 

Models  used  in  our  method  are  identical  to  those  in 
Chin’s  [1].  Following  the  approach  of  the  Deal-Grove  one¬ 
dimensional  oxidation  model  [3],  the  oxidant  diffusion  is 
treated  as  a  steady-state  process  with  uniform  diffusivity. 
The  oxidant  flux  F  is  given  by 

F  =  -DVC  (1) 

where  D  and  C  are  the  diffusion  coefficient  and  concen¬ 
tration  of  oxidant.  The  steady-state  assumption  requires 


V  •  F  =  0  (2) 

this  immediately  suggests  that  C  satisfies  the  Laplace’s 
equation  if  D  is  a  constant; 

V-C  =  0.  (3) 

The  oxide  is  modeled  as  a  viscous  incompressible  fluid 
described  by  the  creeping  flow  and  continuity  equations: 

=  VP  (4) 

V  ■  u  =  0  (5) 


where  ji.  « ,  and  P  are  the  viscosity,  velocity  and  pressure 
of  oxide  respectively.  The  viscous  fluid  approximation  is 
valid  as  long  as  the  oxidation  temperature  is  above  960°C. 
Below  the  critical  temperature,  oxide  exhibits  strong  elas¬ 
tic  characteristics  [4]. 

Incompressible  flow  in  two  dimensions  can  be  described 
by  a  stream  function  which  is  related  to  u  by 

(6) 


dx 


where  .f  and  y  are  unit  vectors  in  the  x  and  _v  directions. 
For  creeping  flow,  satisfies  the  biharmonic  equation: 

V*’!'  =  0  (7) 

where  V*  is  the  Bi-Laplace  operator: 


V  =  ^  + 


-t- 


d.r"  dr-dv-  dv"' 


(8) 


Normally  a  vector,  the  vorticity  w  is  treated  here  as  a  sca¬ 
lar  since  it  is  always  normal  to  the  plane  in  two-dimen¬ 
sional  problems;  it  is  given  by 

dw,.  bu. 
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where  the  subscripts  x  and  y  denote  the  x  and  y  compo¬ 
nents  of  vector  u.  It  follows  from  (6)  that  cj  is  related  to 
*  by 

u  =  (10) 

An  important  feature  in  this  biharmonic  formulation  is  that 
P  and  w  are  related  by  the  Cauchy-Rieman  equations  [5]: 

dcj  _  9(P//l) 

dx  dy 

^ 

dy  dx  '  ^  ^ 

This  conjugate  harmonic  relationship  enables  us  to  obtain 
P  from  'Hr  with  the  following  operator; 

(12) 

dubbed  as  “del-square  conjugate”  operator.  This  operator 
takes  the  Laplacian  of  ^  to  produce  -w,  and  conjugates 
the  result  to  yield  —Pin. 

With  the  principal  components,  namely,  u  and  P,  read¬ 
ily  available  from  all  types  of  boundary  conditions  can 
be  modeled.  A  list  of  those  boundary  conditions  is  given 
in  the  Appendix. 

As  a  masking  material  impermeable  to  most  oxidants, 
silicon  nitride  makes  LOCOS  possible.  In  addition  to 
blocking  the  supply  of  oxidant  species  from  the  ambient, 
it  acts  as  barrier  to  the  flow  of  oxide.  A  simplified  beam 
bending  equation  described  in  [9]  had  been  adopted  to 
model  bending  of  the  silicon  nitride  mask,  but  no  satis¬ 
factory  results  were  obtained.  Thus  for  the  time  being  the 
silicon  nitride  mask  is  assumed  to  be  sufficiently  thin,  that 
is  less  than  300 A,  and  flexible  that  it  does  not  exert  any 
normal  force  on  the  oxide.  It  should  be  pointed  out  that 
different  silicon  nitride  models  can  be  incorporated  into 
our  program  readily  because  they  only  enter  as  boundary 
conditions. 


Fig.  1.  Boundary  segmentation  for  a  SemiRox  structure.  Segments  are  in¬ 
dicated  by  tic  marks.  Three  different  types  of  boundary  are  present:  free 
surface  where  the  oxide  is  in  contact  with  the  gas,  silicon  nitride  mask, 
and  oxide-silicon  interface. 

their  net  potential  contribution  matches  the  prescribed 
boundary  condition,  and  this  gives  the  solution  for  the 
whole  region.  A  detailed  report  of  the  biharmonic  method, 
applicable  to  general  creeping  flow  problems,  will  appear 
in  another  article. 

As  with  simulations  of  other  physical  continua,  we  limit 
our  oxidation  modeling  to  a  small  region  of  a  wafer  due 
to  practical  considerations.  Reflective  boundaries  must  be 
introduced,  one  on  the  left,  one  on  the  right,  to  enclose 
the  region.  However,  artificial  walls  may  not  be  arbitrarily 
used  in  integral  equation  technique,  for  the  sharp  corners 
created  between  the  natural  and  artificial  boundaries  will 
introduce  large  errors  in  the  numerical  solutions.  To  cir¬ 
cumvent  this  problem,  simulation  structures  are  mirror 
imaged  to  become  periodic  symmetric.  In  this  case,  the 
symmetry  lines  represent  reflecting  boundaries.  Dimen¬ 
sions  are  scaled  to  be  Iv  periodic  and  K  is  modified  to  be 
periodic.  As  an  example,  (13)  is  rewritten  as 

C(p)  =  \  a(^  2  K{2-xnx;  p  -  q)  dq 

rt  »  -  00 

/•  oe 

-(-1  o{q)  Yi  Kil-xnx  +  p  -  q)  dq 

JdB.  n  -  -  00 


Solution  Method 

Harmonic  and  biharmonic  problems  can  be  solved  effi¬ 
ciently  using  the  integral  form  of  the  governing  equations. 
Both  the  oxidation  diffusion  and  oxide  flow  problems  are 
modeled  using  the  source  distribution  technique  in  which 
the  solutions  are  expressed  in  terms  of  source  distributions 
on  the  boundaries  [6],  [7]: 

Qp)  =  \  o(q)K(P  -  q)dq  (13) 

JdB 

^(p)  =  I  a^{^K^{p  -  q) 

J3B 

-I-  ai{q)K4p  -  q)  dq  (14) 

where  a’s  and  K's  are  the  sources  and  the  kernels,  re¬ 
spectively.  and  dB  denotes  the  boundaries.  A  simple  anal¬ 
ogy  for  the  harmonic  problem  is  the  potential  distribution 
in  an  enclosed  space.  In  this  case,  o  is  just  positive  and 
negative  charges,  and  K  the  potential  distribution  due  to 
a  point  charge.  Charges  are  arranged  in  such  a  way  that 


o{q)K{p  -  q)  dq 

3Bi 

-F  \  a{q)K(p  -  q)  dq  (15) 

JdB. 

where  dBi  is  a  period  of  the  lower  bundary,  and  of  the 
upper  boundary.  The  tilde  symbol  denotes  periodicity,  that 
is  K(p)  =  K(p  -F  27rat)  for  any  integer  n.  Forms  of  the 
periodic  kernels  are  given  in  the  Appendix. 

The  boundaries  are  divided  into  straight  line  segments, 
as  illustrated  in  Fig.  1.  Typically  there  are  20  to  35  seg¬ 
ments  per  half  period  for  each  of  the  top  and  the  bottom 
boundaries.  The  lower  limit  is  due  to  accuracy  consider¬ 
ation,  and  the  upper  limit  is  constrained  by  computation 
speed.  After  the  boundaries  are  segmented,  the  source 
distributions  on  those  segments  are  approximated  with 
step  functions,  that  is  the  source  strength  is  constant  along 
the  length  of  a  segment.  The  midpoints  of  the  segments 
are  chosen  to  be  the  collocation  points  (test  points)  where 
boundary  conditions  are  imposed.  We  recall  that  bound- 
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ary  parameters  are  obtained  by  applying  differentiation 
operators  on  C  and  'ft.  According  to  (13)  and  (14),  these 
operators  in  turns  operate  on  K,  K^,  and  K^.  Thus  for  the 
oxidant  diffusion  problem,  the  boundary  condition  at 
boundary  segment  i  is  given  by 

N  j. 

bi  =  S  Uj  \  £,(A:(p,  -  qj)]  dqj  (16) 

y«I  J3Bj 

where  i>„  £,,  p,-  are  the  boundary  condition,  associated 
operator,  and  midpoint  of  segment  i,  respectively,  and  N 
is  the  number  of  segments. 

Each  integral  within  the  summation  of  (16)  is  obtained 
from  either  a  fourth-order  or  a  sixth-order  Newton-Cotes 
quadrature  [8] .  The  quadrature  scheme  however  fails  when 
the  source  and  collocation  points  coincide  and  produce  a 
singular  integrand.  In  this  case  the  kernels  are  converted 
into  nonperiodic  form  by  taking  second-order  Taylor  series 
expansions,  as  listed  in  the  Appendix.  The  integral  is  then 
obtained  analytically.  Limits  are  taken  such  that  sources 
always  lie  outside  the  enclosed  region. 

The  system  of  equations  for  all  d,  as  described  by  (16) 
can  be  arranged  into  a  matrix  format, 

Aa  =  b  (17) 

from  which  the  values  for  <r’s— the  unknowns— are  ob¬ 
tained  with  a  suitable  technique  such  as  Gaussian  elimi¬ 
nation.  The  size  of  <4  is  N  by  N  for  the  oxidant  diffusion 
and  2N  by  2N  for  the  oxide  flow.  Equation  (16)  is  used 
again,  this  time  with  known  a’s  and  different  operators, 
to  obtain  boundary  parameters  desired.  Note  that  we  are 
free  to  choose  any  operators  as  long  as  they  are  reason¬ 
able.  Our  technique  thus  lends  itself  to  extracting  bound¬ 
ary  parameters  such  as  oxidant  flux  and  shear  stress  at  the 
oxide-silicon  interface. 

The  freedom  to  define  operators  for  extracting  param¬ 
eters  from  the  intermediary  solution  is  a  unique  feature  of 
our  approach.  In  contrast,  the  choice  of  parameters  that 
can  be  obtained  from  Green’s  third  identity  (Green’s  for¬ 
mula),  which  is  used  in  [1]  and  [2],  is  very  limited. 
Granted  Green’s  third  identity  does  not  require  interme¬ 
diary  solutions,  that  is  the  sources,  and  thus  is  probably 
more  accurate;  ours,  on  the  other  hand,  is  more  flexible. 

A  first-order  method  is  used  in  the  time-stepping 
scheme.  After  growth  rate  is  computed  from  the  solution 
for  the  oxidant  diffusion,  the  oxide  flow  problem  is  solved 
to  determine  the  velocities  of  the  boundary  segments.  The 
boundary  segments  are  then  moved  according  to  the  time- 
step  size  chosen.  This  cycle  of  operations  is  repeated. 

Applications 

We  now  discuss  a  few  examples  on  how  to  use  our 
method.  These  technological  applications  include  simula¬ 
tion  of  LOCOS  processes,  analysis  of  stress  associated 
with  deformation,  and  optimization  of  ambient  pressure 
to  reduce  oxidation  time. 


Fig.  2.  Simulation  of  a  SemiRox  structure  with  pad  oxide  (relief  oxide)  of 
500  A .  Final  field  oxide  thickness  is  approximately  SOOO  A .  (a)  Outline 
of  the  oxide  layer  before  and  after  oxidation,  (b)  Normal  stress  distri¬ 
bution  at  the  oxide-silicon  interface.  Note  that  the  positive  direction  is 
into  the  oxide  bulk,  (c)  Shear  stress  distribution  at  the  oxide-silicon  in¬ 
terface.  Direction  of  force  is  tangent  to  the  interface. 


A.  LOCOS  Simulations 

The  two  LOCOS  presented  are  the  semi-recess  oxida¬ 
tion  (SemiRox)  and  sidewall  masked  isolation  (SWAMI) 
techniques. 

The  SemiRox  structure  has  a  pad  oxide  of  500  A  and  a 
thin  silicon  nitride  mask.  As  discussed  earlier,  the  mask 
is  assumed  to  be  sufficiently  thin  so  that  is  does  not  gen¬ 
erate  significant  forces  when  bent.  Oxidation  is  carried 
out  at  1025 °C  in  wet  Oi.  The  shape  of  the  oxide  layer  is 
shown  in  Fig.  2(a)  with  normal  and  tangential  components 
of  the  stress  distribution  along  the  oxide-silicon  interface 
given  in  Fig.  2(b)  and  (c),  respectively.  The  value  for  the 
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Fig.  3.  Sloped-wall  SWAMI  oxidation.  Pad  oxide  is  200  A. 

oxide  viscosity  is  taken  from  [9]: 

M7)  =  /*o  exp  07) 

where  fio  =  1.586  X  10"'°  poises,  =  5.761  eV,  k  is  the 
Boltzmann  constant,  and  Tis  the  temperature.  The  normal 
stress  distribution  is  dipole  due  to  the  highly  nonuniform 
oxide  growth  rate  at  the  edge  of  the  silicon  nitride  mask. 
Whether  the  stress  is  compressive  (negative)  or  tensile 
(positive)  depends  on  the  oxide  is  being  slow  down  or 
pulled.  As  the  oxide  layer  gets  thicker,  the  growth  rate 
decreases  and  becomes  more  uniform,  hence  the  decrease 
in  the  normal  stress.  Being  more  complicated,  the  behav¬ 
ior  of  the  shear  stress  cannot  be  described  in  intuitive 
terms. 

For  the  SWAMI  process,  we  simulate  an  improved  tech¬ 
nique  that  uses  sloped  wall  [10].  Simulation  parameters 
are  200A  pad  oxide,  975‘’C,  and  wet  Oj.  The  sidewall  is 
sloped  at  54.74®;  silicon  nitride  covers  the  top  and  side- 
wall,  and  extends  slightly  into  the  trench.  Fig.  3  shows 
the  initial  and  Anal  shapes  of  the  structure. 

In  general,  our  results  from  different  simulations  agree 
well  Chin’s.  The  discrepancy  (approximately  20  to  30 
percent)  in  the  normal  stress  is  mainly  due  to  the  different 
ways  of  calculating  stress— Chin  approximated  stress  with 
hydrodynamic  pressure  whereas  we  dealt  with  the  full 
stress  tensor. 

B.  Deformation  Stress  Calculation 

The  plastic  yield  strength  of  single-crystal  silicon  is  a 
strong  function  of  temperature.  At  room  temperatures, 
single-crystal  silicon  is  known  to  fracture  but  never  de¬ 
form  when  subjected  to  stress  [11].  Enhanced  oxidation 
notwithstanding,  defects  do  not  appear  in  plasma-assisted 
oxidation,  which  is  carried  out  at  moderately  elevated 
temperatures  (between  300°C  and  600“C)  [12].  It  is 
mainly  in  conventional  high-temperature  thermal  oxida¬ 
tion  that  defects  are  observed. 

Slip  is  the  major  defect  species,  its  occurrence  is  con¬ 
fined  to  the  {111 )  planes.  On  a  { 100}  wafer  typically  used 
in  making  MOS  devices,  the  {111}  planes  intersect  on  the 
surface  at  right  angle  relative  to  one  another.  Devices  are 
usually  aligned  with  these  “fault  lines’’  for  ease  of  dicing. 


Fig.  3.  Shear  stress  in  slip  planes.  The  shape  of  the  oxide  layer  is  shown 
in  (he  insert,  together  with  the  slip  directions.  Simulation  structure  is  the 
sloped-wall  SWAMI  oxidation  given  in  Fig.  3. 

If  we  look  at  a  suitable  cross  section  of  a  device,  two  { 111 } 
planes  make  54.74®  with  the  surface,  as  shown  in  Fig.  4. 
Slip  occurs  when  the  shear  stress  on  these  slip  planes  ex¬ 
ceeds  the  plastic  yield  strength  of  silicon. 

The  advantage  of  source  distribution  technique  is  ap¬ 
parent  now.  With  the  solutions  defined  by  the  intermedi¬ 
ary  sources,  the  normal  stress  and  the  shear  stress  tangent 
to  the  oxide-silicon  interface  can  be  determined.  The  same 
sources  can  also  be  used  to  calculate  the  shear  stress  in 
the  slip  planes.  This  is  done  by  redefining  the  local  coor¬ 
dinate  ij  —  5  of  the  shear  stress  operator  to  align  with  the 
slip  planes.  Fig.  5  shows  the  deformation  stress  associated 
with  the  two  slip  planes  for  the  SWAMI  structure  men¬ 
tioned  easier. 

Strictly  speaking,  we  cannot  accurately  model  defor¬ 
mation  shear  stress  in  two  dimensions;  the  reason  has  to 
do  with  the  slip  directions.  Consider  one  of  the  { 111 }  plane 
shown  in  Fig.  4.- There  are  three  <110)  slip  directions  with 
this  slip  plane,  one  normal  to  the  paper,  the  other  two  in¬ 
tersecting  the  paper  at  a  low  angle  (30®).  Our  method  can 
only  approximate  the  shear  stress  associated  with  the  lat¬ 
ter  two,  while  ignoring  that  of  the  first  one. 

C.  Ramped-Pressure  Oxidation 
In  the  third  application  the  feasibility  of  reducing  oxi¬ 
dation  time  by  ramping  the  ambient  pressure  is  investi¬ 
gated.  To  produce  the  required  oxide  thickness,  the  field 


19Sg 


IEEE  transactions  ON  ELECTRON  DEVICES.  VOL.  ED-32.  NO.  10.  OCTOBER  1983 


Conclusion 

A  boundary  integral  equation  method  has  been  devel¬ 
oped  to  model  oxidation  diffusion  in  two  dimensions.  It  is 
a  simple  technique  due  to  the  fact  that  the  oxidation  dif¬ 
fusion  can  be  modeled  as  a  harmonic  problem  and  the  ox¬ 
ide  motion  as  a  biharmonic  problem.  In  theory,  problems 
such  as  the  effect  of  stress  on  surface  reaction  rate  can  be 
readily  incorporated  into  our  method  as  long  as  they  can 
be  treated  as  a  change  in  the  boundary  conditions. 

OXIDATION  TIME  (hr*) 


3  ' 

3 

3 


1 


-  RAMKO  PRESSURE 

- ATMOSPHERIC  PRESSURE 


(a) 


OXIDATION  TIME  (hr*) 


(b) 

Fig.  6.  (a)  Peak  stress  distribution  (in  slip  planes)  of  conventional  and 
ramped-pressure  oxidation,  plotted  as  a  function  of  time.  Simulation 
structure  is  similar  to  the  SemiRox  oxidation  shown  in  Fig.  2.  (b)  Am¬ 
bient  pressure  curve  for  ramped-pressure  oxidation. 


oxidation  step  must  subject  a  wafer  to  high  temperatures 
for  long  periods,  thus  causing  significant  impurity  redis¬ 
tribution  in  the  substrate.  Knowing  that  stress  is  propor¬ 
tional  to  oxidation  rate  and  that  oxidation  rate  decreases 
as  oxide  layer  gets  thicker,  we  can  attempt  to  shorten  the 
oxidation  time  by  gradually  ramping  the  ambient  pressure 
beyond  atmospheric.  In  doing  so  we  may  however  gener¬ 
ate  a  stress  that  exceeds  the  plastic  yield  strength  of  sili¬ 
con.  hence  damaging  the  substrate.  The  plastic  deforma¬ 
tion  can  be  avoided  if  we  increase  the  ambient  pressure  in 
such  a  way  that  the  maximum  stress  always  stays  below 
the  critical  limit. 

The  test  structure  used  in  this  study  is  SemiRox,  with 
pad  oxide  of  500 A.  Oxidation  is  carried  out  at  950'’C  in 
wet  O2;  the  final  field  oxide  thickness  is  approximately 
8000  A .  A  simulation  is  first  done  at  constant  atmospheric 
pressure  to  decide,  in  an  arbitrary  manner,  a  safe  value 
for  the  maximum  shear  stress  that  the  silicon  substrate  may 
withstand  in  its  slip  planes.  This  value  is  then  used  as  a 
ceiling  in  a  second  simulation  which  ramps  the  ambient 
pressure.  Our  result  indicates  that  the  oxidation  time  can 
be  reduced  by  half,  roughly.  Fig.  6(a)  shows  the  maximum 
shear  stress  as  a  function  of  time  for  the  case  when  am¬ 
bient  pressure  is  maintained  at  atmospheric  and  for  the 
case  when  it  is  ramped.  A  plot  of  the  ambient  pressure  in 
the  ramped-pressure  oxidation  is  given  in  Fig.  6(b),  the 
pressure  begins  at  1  atm  and  increases  to  4  towards  the 
end. 


Appendix 

Boundary  Parameters  and  Operators  for 
Creeping  Flow 


Boundary  parameters  that  may  be  obtained  from  the 
stream  function  are  as  follows: 

(1)  Normal  component  of  flow  velocity: 


u 


(Al) 


(2)  Tangential  component  of  flow  velocity: 


(A2) 

=  V 

(3)  Vorticity: 

u  = 

(A3) 

(4)  Hydrodynamic  pressure: 

P  = 

(A4) 

(5)  Normal  stress: 

=  ^  V^- 


dvdl 


(6)  Shear  stress: 


-  Ml 

\d^  dij/ 

■  a^  a^  ■ 

=  M 

.a'r? 

'ir 


(A5) 


(A6) 


where  ij  is  unit  normal  vector  pointing  into  the  oxide  bulk 
and  I  is  unit  tangent  vector  90“  clockwise  from  ij. 


Boundary  Conditions  for  LOCOS 

At  the  free  surface  the  oxidant  concentration  is  main¬ 
tained  at  the  equilibrium  concentration  C*;  both  the  nor¬ 
mal  and  shear  stress  components  are  zero  since  there  is 
no  force  exchange  between  the  oxide  and  the  gas: 

C  =  C* 

T,(  =  0.  (A7) 
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At  the  oxide-silicon  interface  the  oxidation  rate  is  pro¬ 
portional  to  both  the  oxidant  concentration  and  nor- 
marcomponent  of  oxidant  flux,  thus  the  following  form  of 
boundary  condition  is  applied; 
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D  —  =  k,C  (ASA) 

dry 

where  k,  is  the  chemical  surface  reaction  rate.  The  oxide- 
silicon  interface  recesses  into  the  substrate  at  a  velocity 

ksC  . 

H  r,  (ASB) 

while  oxide  is  injected  into  the  bulk  at 

u  =  {a  -  l)~fi.  (A8C) 

ajVi 

o  is  the  volume  of  oxide  produced  from  a  unit  volume  of 
silicon,  it  is  approximately  2.24.  is  the  number  of  ox¬ 
idant  molecules  required  to  form  a  unit  volume  of  oxide. 

At  the  oxide-silicon-nitride  interface,  the  oxidant  flux 
has  a  zero  normal  component; 

7-  =  0.  (A9A) 

di} 

The  silicon  nitride  mask  is  assumed  to  be  sufficiently  thin 
such  that  is  does  not  exert  normal  stress  on  the  oxide; 

r„  =  0  (A9B) 

while  maintaining  a  no-slip  boundary  condition; 

uj  =  0.  (A9C) 

Kernels 

Among  the  3  periodic  kernels,  only  cannot  be  ex¬ 
pressed  in  simple  functions,  but  its  derivatives  can. 

=  J  log  [2(cosh  (y)  -  cos  (.x))] 

=  ^v  log  [2(cosh  (y)  -  cos  (.r))] 

=  yK(r) 

dy 

=  y^^(r).  (AlO) 

The  corresponding  nonperiodic  Taylor  series  approxima¬ 
tions  are  obtained  by  expanding  the  trigonometric  and  hy¬ 
perbolic  functions  to  second  order  at  r  =  0. 

K(r)  =  i  log  (.r  +  jr) 

=  i  y  log  (jr  +  y) 

Kz(r)  =  i  X  log  (jr  +  y)  -  X.  (All) 


K(r) 

^.(r) 

f  ^2(^0 

dx 

dy 
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ABSTRACT 


The  scheduling  of  a  VLSI  manufacturing  facility  is  being 
studied  using  MIT’s  new  Integrated  Circuits  Laboratory  (ICL)  as  a 
testbed.  This  paper  describes  the  facility,  outlines  the  schedu¬ 
ling  methodology  that  is  under  investigation  and  presents  the 
design  of  the  initial  preliminary  scheduling  system. 

INTRODUCTION 

Scheduling  is  the  process  of  allocating  the  use  of  a  set  of 
scarce  resources  over  time  so  as  to  best  meet  a  set  of  objec¬ 
tives.  In  this  case,  the  resources  are  the  machines,  techni¬ 
cians,  and  support  systems  of  a  semiconductor  fabrication  facili¬ 
ty.  The  objectives  are  to  satisfy  the  needs  of  the  users  of  the 
facility  with  minimal  delay.  Because  the  resources  are  finite, 
there  are  times  when  users’  needs  may  conflict.  There  must  be  a 
set  of  policies  to  deal  with  these  conflicts,  and  the  scheduler 
must  resolve  them  in  a  way  that  is  consistent  with  system  policy. 

THE  MIT  INTEGRATED  CIRCUITS  LABORATORY 

The  MIT  Integrated  Circuits  Laboratory  (ICL)  is  a  new  faci¬ 
lity  for  semiconductor  research.  It  will  also  serve  as  both  a 
customer  and  laboratory  for  scheduling  algorithm  development. 

The  building  contains  offices,  clean  rooms,  laboratory  space,  and 
maintenance  and  support  areas.  Clean  rooms  have  low  particle 
counts  due  to  constant  flow  of  filtered  air.  Clean  rooms  contain 
either  communal  or  private  equipment.  Laboratories  are  not  ne¬ 
cessarily  located  in  clean  rooms. 

House  systems  are  threaded  through  the  building  providing 
services  such  as  air  flow,  humidity  control,  gaseous  nitrogen, 
deionized  water,  and  solvent  waste  removal. 

A  network  of  computers  and  terminals  provides  the  hardware 
needed  by  the  Computer-Aided  Fabrication  (CAF)  project.  A  goal 
of  this  project  is  the  efficient  use  of  computers  in  semiconduc¬ 
tor  fabrication.  As  a  part  of  the  CAF  project,  scheduling  is 
under  study  for  both  coordinating  the  activities  in  the  ICL  (to 
make  the  best  possible  use  of  valuable  resources  such  as  ma¬ 
chines,  technicians,  power,  steam,  chilled  water,  and  chemicals) 
and  to  develop  algorithms  for  industry.  The  purpose  of  this 
paper  is  to  describe  our  preliminary  work  in  scheduling  the  ICL. 

The  ICL  facility  is  run  by  a  team  of  professionals.  Techni¬ 
cians  maintain  and  operate  equipment.  Process  engineers  are 
responsible  for  groups  of  machines  and  the  processes  that  run  on 
those  machines.  A  manager  is  responsible  for  the  smooth  opera¬ 
tion  of  the  facility. 

The  facility  is  used  by  about  100  students  and  20  faculty 
members  whose  interests  include  submicron  structures,  semiconduc¬ 
tor  materials,  semiconductor  devices,  semiconductor  processes, 
and  CAF  research. 


- 
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Most  research  involves  perforning  operations  on  silicon 
wafers.  A  group  of  wafers  that  is  processed  together  is  called  a 
lot.  The  sequence  of  activities  required  to  produce  a  lot  is 
called  a  flow.  A  linear  sequence  of  operations  is  called  a 
liS§§£  A  flow  can  branch.  Lots  can  be  split  for  rework  or 

experiaentation.  The  flow  of  a  research  lot  may  also  depend  on 
■easureaents  aade  after  a  process  step.  Wafers  are  discarded  if 
unusable. 

In  addition  to  production  and  laboratory  operations,  other 
activities  such  as  maintenance,  training,  inspection,  repair,  and 
aachine  setups  must  be  scheduled. 

PROBLEM  DESCRIPTION 

The  scheduling  problem  is  to  choose  a  vector  T  where  TJ  is 
the  time  that  activity  number  J  is  scheduled  to  begin.  Tj  is 
selected  from  the  set  of  legal  times  for  activity  J.  This  set  is 
limited  due  to  personal  schedule  constraints,  aachine  repair 
states,  and  previously  scheduled  activities  such  as  preventive 
maintenance.  T  is  further  constrained  by  limited  capacities  of 
machines,  house  systems,  and  technicians,  and  by  radio  interfe¬ 
rence  among  machines  and  by  conflicting  setup  requirements. 

METHODOLOGY 

The  methodology  selected  for  studying  the  scheduling  of  this 
facility  involves  the  following: 

1)  Study  semiconductor  fabrication  technology  and  the  ICL  to 
determine  important  features  of  the  problem.  The  schedule  is 
constrained  by  the  structure  of  the  system.  For  example: 

a)  Exposing  wafers  to  air  can  dramatically  reduce  yield. 
Therefore  certain  sequences  of  activities  must  not  have  any 
delay . 

b)  The  house  system  which  supplies  deionized  water  has  a  300 

gallon  storage  capacity,  and  a  30  gallon  per  minute  generation 
capacity.  A  single  wet  sink  uses  20  gallons  per  minute. 

An  effective  schedule  must  limit  use  of  deionized  water  to  within 
capacity. 

c)  The  ion  implanter  has  two  different  types  of  setups.  The 
source  and  the  wafer  holder  can  both  be  changed.  The  scheduling 
system  must  insure  that  both  setups  are  appropriate,  or  it  must 
reserve  machine  and  technician  time  to  change  them. 

2)  Specify  objectives  of  scheduling 

Possible  scheduling  goals  include  minimizing  cycle  time  (the 
time  a  wafer  spends  in  the  system),  maximizing  throughput,  mini¬ 
mizing  cost,  or  satisfying  a  preferred  customer.  Low  cycle  time 
requires  small  buffers,  high  throughput  requires  larger  buffers. 


Lower  costs  means  things  should  wait  until  several  batches  can  be 
processed  together.  Satisfying  a  preferred  customer  means  cer¬ 
tain  lots  move  first  no  matter  what  the  consequences.  Bach  of 
these  objectives  can  lead  to  a  different  schedule. 

The  initial  project  goal  is  to  provide  a  simple,  effective 
scheduling  system  for  the  MIT  ICL  Later  versions  of  the  sche¬ 
duling  algorithm  will  meet  criteria  such  as  those  described 
above . 

To  use  optimal  control  techniques  to  help  develop  feedback 
rules  for  scheduling  requires  a  numerical  objective  function  L. 
One  possible  objective  is  the  minimization  of  the  average  expec¬ 
ted  time  between  the  arrival  of  an  order  for  a  part  and  the 
completion  of  the  part. 

3)  Develop  a  solution  technique. 

The  scheduling  approach  will  follow  the  framework  described 
in  Gershwin  et  al.  (1986).  In  that  framework,  the  capacity 
constraints  of  the  manufacturing  system  are  explicitly  modeled 
and  feedback  control  strategies  are  developed  to  respond  to 
random,  potentially  disruptive  events 

Because  of  the  complexity  of  the  problem,  hierarchical  sche¬ 
duling  methods  are  being  developed.  The  top  level  of  the  hierar¬ 
chy  deals  with  relatively  long  term  issues  such  as  maintenance 
frequencies,  setup  frequencies,  historical  demand,  etc.  It  allo¬ 
cates  time  on  each  of  the  machines  for  each  class  of  activities. 
Since  there  is  a  substantial  setup  time  between  certain  activi¬ 
ties  (such  as  when  there  are  two  different  kinds  of  gases  or 
impurities  in  a  chamber) ,  the  frequencies  of  the  changeovers  must 
be  chosen  with  care. 

For  example,  given  the  best  current  estimates  about  usage 
during  a  given  future  time  period  (such  as  a  week  three  weeks 
hence)  the  top  level  of  the  scheduler  may  decide  that  the  ion 
implanter  will  be  set  up  with  the  hearth  stage  from  9:00  AM 
Monday  until  11:00  AM  Wednesday,  after  which  it  will  be  changed 
over  to  the  cold  stage.  Maintenance  may  then  be  scheduled  until 
2:00  PM,  followed  by  a  setup  with  the  hearth  stage.  These  allo¬ 
cations  are  chosen  on  the  basis  of  historical  data,  without 
reference  to  current  requests  for  that  machine  during  that  time 
period.  In  addition,  some  amount  of  time  is  built  into  the  setup 
intervals  to  allow  for  such  disruptions  as  machine  failures. 

As  actual  demand  data  arrives,  the  lowest  level  of  the  hier¬ 
archical  scheduler  assigns  requests  to  these  time  intervals.  For 
example,  the  period  from  9:00  Monday  until  11:00  Wednesday  of 
that  week  is  only  used  to  satisfy  requests  for  the  ion  implanter 
with  the  hearth  stage.  If  the  demand  for  the  ion  implanter  with 
the  hearth  stage  is  greater  than  the  allocated  time,  or  if  some 
users  have  schedule  conflicts  that  prevents  their  use  during  that 
period,  the  excess  is  moved  to  the  interval  starting  at  2:00 
Wednesday . 


The  Biddle  level  of  the  scheduler  has  two  functions:  first, 
it  responds  to  random  events  such  as  aachine  failures.  When  a 
machine  goes  down,  it  adjusts  the  schedule  in  accordance  with 
facility  policies  about  fairly  allocating  the  inconvenience 
caused  by  the  failure  among  the  users  that  are  affected.  Second, 
it  adjusts  the  setup  intervals.  If,  for  example,  during  some 
week  there  is  greater  than  usual  demand  with  the  hearth  stage  and 
less  for  the  cold  stage,  the  11:00  Wednesday  changeover  may  be 
deferred  until  noon. 

If  there  is  a  systematic  discrepancy  between  the  presche¬ 
duled  setup  intervals  and  the  demand,  it  is  the  responsibility  of 
the  highest  level  of  the  algorithm  to  adjust  the  lengths  of  the 
intervals . 

This  preliminary  statement  is  far  from  complete.  We  antici¬ 
pate  using  methods  such  as  those  of  Kimemia  and  Gershwin  (1983) 
and  Gershwin,  Akella,  and  Choong  (1985)  to  respond  to  machine 
failures,  and  Graves  et  al.  (1983)  and  Watts  (1986)  for  the 
detailed  dispatching  of  lots. 

4)  Develop  an  appropriate  database. 

The  detailed  schedule  is  strongly  dependent  on  a  wide  varie¬ 
ty  of  kinds  of  information.  Required  data  structures  must  repre¬ 
sent  flows,  personal  schedules,  machine  schedules,  lot  schedules, 
reservations,  machine  maintenance,  machine  setup  states,  machine 
repair  states,  the  effects  of  failures  of  specific  house  systems, 
the  maintenance  of  house  systems,  the  priorities  of  lots,  the 
facilities  authorization  list,  and  the  skills  list. 

5)  Teat  the  entire  system 

The  scheduling  algorithms  developed  will  be  used  to  run  the 
ICL  as  an  experimental  test  of  these  ideas. 

VERSION  1  SCHEDULING  SYSTEM 

As  a  first  step  towards  solving  the  problem, 
duling  system  is  being  implemented.  This  system 
be  little  more  than  a  computerized  sign-up  sheet, 
are  to 

a)  Allow  users  to  reserve  equipment. 

b)  Inform  users  when  machines  are  down,  scheduled  for  mainte¬ 
nance,  or  reserved  by  someone  else. 

The  table  containing  the  desired  information  is  shown  in  Fig 
1.  There  is  one  column  for  every  resource.  The  rows  correspobd 
to  time  slots.  Future  schedules  are  listed  below  the  current 
line,  previous  ones  above.  The  user  sees  a  window  into  this 
table  as  shown  in  Fig.  2.  Because  of  limited  screen  size,  only 
four  columns,  representing  four  resources,  are  shown  at  a  time. 


a  simple  sche- 
is  intended  to 
Its  objectives 


CONCLUSION 


Scheduling  VLSI  facilities  is  a  new  and  important  field. 
Competition  is  pushing  American  manufacturers  away  from  mass 
production  and  into  Application-Specific  Integrated  Circuits 
(A.S.I.C.).  The  rising  coat  of  new  facilities  is  putting  pres¬ 
sure  to  process  research  lots  on  the  same  line  as  production 
lots.  Both  of  these  trends  increase  the  number  of  flows  in  a 
single  facility.  More  flows  make  scheduling  more  difficult  and 
thus  increase  the  opportunity  for  profiting  from  scheduling. 

This  will  make  VLSI  scheduling  an  important  activity  in  the  near 
future . 
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